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Abstract 

The analytical solutions of the equations describing impurity diffusion due to migra- 
tion of nonequilibrium impurity interstitials were obtained for the impurity redistribution 
during ion implantation at elevated temperatures and for diffusion from a doped epitax- 
ial layer. The reflecting boundary condition at the surface of a semiconductor and the 
conditions of constant concentrations at the surface and in the bulk of it were used in 
the first and second cases, respectively. On the basis of these solutions hydrogen diffusion 
in silicon during high-fluence low-energy deuterium implantation and beryllium diffusion 
from a doped epi-layer during rapid thermal annealing of InP/InGaAs heterostructures 
were investigated. The calculated impurity concentration profiles agree well with exper- 
imental data. The fitting to the experimental profiles allowed us to derive the values of 
the parameters that describe interstitial impurity diffusion. 
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1 Introduction 

The recent years numerical methods have been widely used for simulation of solid state 
diffusion of ion-implanted dopants (see, for example, [H [2]). As a rule, to simulate the 
impurity diffusion a system of equations describing a coupled diffusion of different mo- 
bile species and their quasichemical reactions during annealing is solved. Due to a great 
number of differential equations and the complexity of the system as a whole, the prob- 
lem of the correctness of a numerical solution is very important. One of the best ways 
to verify the correctness of an approximate numerical solution is a comparison with the 
exact analytical solution of the boundary value problem under consideration. Such an- 
alytical solutions can be derived for the special simplest cases of dopant or point defect 
diffusion processes. For example, in Ref. [3] an analytical solution for the point defect 
diffusion based on the method of Green's functions was obtained. It was supposed in [3] 
that nonequilibrium point defects were continuously generated during ion implantation of 
impurity atoms and diffused to the surface and into the bulk of a semiconductor. The sur- 
face was considered to be a perfect sink for point defects. In Ref. [1] a process of impurity 
diffusion during ion implantation at elevated temperatures was investigated analytically. 
It was supposed that the implantation temperature was too low to provide a traditional 
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diffusion by the "dopant atom - point defect" pairs, but was enough for the diffusion 
of nonequilibrium interstitial impurity atoms to occur. Unhke [5], in Ref. [3] a system 
of equations, namely, the conservation law for substitutionally dissolved impurity atoms 
and equation of diffusion-recombination of nonequilibrium interstitial impurity atoms 
have been solved analytically by the method of Green's functions. Reflecting boundary 
condition at the surface of a semiconductor has been chosen to describe the interaction of 
interstitial impurity atoms with the interface. Due to this condition, a diffusion problem 
has become symmetric with respect to the point x = 0. For simplicity, the condition of 
zero impurity concentration for x ±oo has been used. It is interesting to note that 
analytical solutions for gold diffusion in silicon due to Frank- Turnbull and due to kick-out 
mechanisms were obtained in Ref. |5] and Refs. O [6], respectively. It was supposed 
that there was a local equilibrium between substitutionally dissolved impurity atoms, va- 
cancies (or self-interstitials for the kick-out mechanism) and interstitial impurity atoms. 
The case of nonequilibrium interstitial impurity atoms was not considered in these pa- 
pers. The very interesting case of coupled diffusion of vacancies and self-interstitials was 
investigated in [TJ [8]. The equations of the diffusion of vacancies or of self-interstitials 
are similar to the equation of diffusion of impurity interstitials. However, the solutions 
obtained in [71 E] are difficult to use for describing the impurity diffusion governed by 
nonequilibrium interstitial impurity atoms, because a condition of local equilibrium was 
used in all these papers. Besides, a generation rate was assumed to be equal to zero in [8J 
or equal to constant value in [7]. Thus, it is reasonable to derive an analytical solution 
for the case of diffusion due to migration of nonequilibrium impurity interstitials. 

The main goal of this paper is to continue the investigation [3] to obtain other ana- 
lytical solutions and compare these solutions with the experimental data. 

2 Original equations 

It is supposed that the processing temperature is too low to provide a diffusion of substitu- 
tionally dissolved impurity atoms, but is enough for the diffusion of impurity interstitials. 
The generation of nonequilibrium interstitial impurity atoms can occur due to ion im- 
plantation, including the case of ion implantation at elevated temperatures, or due to the 
replacement of the impurity atom by self-interstitial from the substitutional position to 
the interstitial one (Watkins effect [9J), or due to dissolution of the clusters that incorpo- 
rate impurity atoms etc. It is also supposed that the impurity concentration in the doped 
regions formed due to migration of nonequilibrium impurity interstitials is smaller than 
77-2 or approximately equal to rij or that impurity atoms in interstitial position are neutral. 
Here Uj is the intrinsic carrier concentration at the processing temperature. Then, the 
system of equations describing the evolution of impurity concentration profiles includes 



(i) a conservation law for substitutionally dissolved impurity atoms: 





(ii) an equation of diffusion for nonequilibrium interstitial impurity atoms: 
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Here C and C^^ are the concentrations of substitutionally dissolved impurity atoms 
and nonequihbrium impurity interstitials, respectively; G^^ is the rate of introducing of 
impurity atoms, which immediately occupy the substitutional positions; d^^ and t^^ are 
the diffusivity and average lifetime of nonequihbrium interstitial impurity atoms, respec- 
tively; G^^ is the generation rate of interstitial impurity atoms. We use a steady-state 
diffusion equation for impurity interstitials, because of the large average migration length 
of nonequihbrium interstitial impurity atoms {Iai ^ Ifaih where Ifaii is the characteristic 
length of the decrease in the impurity concentration) and due to the small average lifetime 
of nonequihbrium impurity interstitials {tai tp, where is the duration of thermal 
treatment). 

The system ([1]), ([2]) or ([1]), ([3]) describes impurity diffusion due to migration of 
nonequihbrium interstitial impurity atoms. To solve this system of equations, appro- 
priate boundary conditions are need. Let us consider, in contrast to [1], the finite-length 
one-dimensional (ID) domain [0,xb], i.e., the domain used in ID numerical modeling, 
and add the following boundary conditions to Eq. ([3]): 
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(6) 



as well as the initial conditions: 



C(x, 0) = Co(x) , G^\x, 0) = C,^/ = const . (7) 

Here, wf , wf , wf and wf , wf , wf are constant coefficients; G^^ is the equilibrium 
value of concentration of interstitial impurity atoms (it is supposed that C^J is equal to 
zero for many cases under consideration). 

To derive an analytical solution of this boundary value problem, the method of Green's 
function [lOj can be used. 



3 Analytical method and solutions 

The suggestion about the immobile substitutionally dissolved impurity atoms allows one 
to solve Eq. ([1]) independently of Eq. ([2]) or Eq. ([3]): 
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We will supplement expression (|8]) with a steady-state solution of Eq. ([3]) obtained by 
the method of Green's functions llOi: 
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Here is the Green's function for Eq. ([3]). Using the standardizing function 

w{x,t) [lO] allows one to reduce the previous boundary value problem to the boundary 
value problem with zero boundary conditions: 
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The Green's function for Eq. ^ with boundary conditions flTTl) and f[T^ has the 
following form [TO] : 
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K = Q[ (x) Q2 (x) - Qi (x) Q2 (a;) = const . (14) 
Here Qi and (^2 are the linearly independent solutions of the homogeneous equation 
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with the following initial conditions on the left boundary: 



Ql(0) = wf(i^^ g'^(0) = -w 



and on the right one: 
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Following [To], one can write the functions ws{x) and wb{x) as 
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Thus, for the reflecting boundary condition 

wf = 1, wf = 0, wf = 
at the surface of a semiconductor {x — 0) and for Dirichlet boundary condition 
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in the bulk {x = xs) the solutions Qi and Q2 have the form 
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Let us consider the process of ion implantation in a semiconductor at an elevated 
temperature (the so-called "hot" implantation). It is established experimentally that 
the main part of the implanted impurity atoms occupies substitutional positions near 
the places where they were stopped. Let us suppose that the remaining nonequilibrium 
atoms occupy interstitial sites and can diffuse during implantation before they transfer to 
the substitutional position. To describe the spatial distributions of both impurity atoms 
directly occupying substitutional position and impurity interstitials generated during im- 
plantation, Gaussian distributions can be used: 
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where is the number of impurity atoms being introduced per unit area per second by 
ion implantation; p^^ is a part of the implanted impurity atoms occupying interstitial 
sites; Rp and ARp are the average projective range of implanted ions and straggling of 
projective range, respectively. 

Taking into consideration expressions (125!) and (!27|) yields 
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Substituting Green's function (j24l) and standardizing function (!28ll into expression ([9]) 
allow one to obtain a spatial distribution of diffusing interstitial impurity atoms: 
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Calculating the integrals on the right-hand side of expression (1291) one can obtain an 
explicit expression for the distribution of interstitial impurity atoms: 
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In a similar way, one can obtain a solution for the case of impurity interstitial recom- 
bination at the surface of a semiconductor. Let us consider, for example, a buried layer 
uniformly doped by impurity atoms. If the impurity concentration is high, generation 
of nonequilibrium interstitial impurity atoms is possible within this layer during thermal 
treatment. The boundary conditions assuming the constant impurity interstitial concen- 
trations at the surface and in the bulk of a semiconductor (Dirichlet boundary conditions) 
can be enforced to describe the interstitial migration 
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where Cg^ is the concentration of interstitial impurity atoms at the surface (it is quite 
likely that C^^ = 0). 

Let us suppose that the generation of impurity interstitials is described by the following 
function: 
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where Qm = const and x^ and xr are the left and right boundaries of the doped layer, 
respectively. 

Then, the Green's function and standardizing function are respectively 
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Substituting pil) . fH2|) . and (143]) into expression ([9]), we obtain a spatial distribution 
of the diffusing interstitial impurity atoms: 
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To calculate the integrals in expression (H5l) . let us consider the following three 
i) If X < xl, then 
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Expressions (H7I) . (HSj) . and (H^ are the partial steady-state solution of the boundary- 
value problem under consideration for the case of step distribution of interstitial generation 
rate fl41l) and zero concentrations of nonequilibrium impurity interstitials on the left and 
right boundaries of the solution domain. Combining these expressions with solution fH6|) 
of the appropriate homogeneous diffusion equation, we obtain the distribution of impurity 
interstitial concentration for the case of the arbitrary boundary concentrations of impurity 
interstitial atoms (see expression 



4 Simulation 

Using the system of Eqs. ([T]) and ^ and appropriate analytical solutions, one can ver- 
ify the correctness of numerical calculations. This system can be also used for modeling 
hydrogen diffusion in silicon. For the case of hydrogen diffusion, the quantities C and 
C^^ are the concentration of trapped (immobile) hydrogen atoms and concentration of 
the neutral mobile hydrogen interstitials, respectively. As example, in Fig. 1 the results 
of simulation of hydrogen diffusion in silicon obtained on the basis of analytical solution 
dHD and (15(1 are presented. For comparison, the experimental data of [TT] are used. In 
[TT] deuterium was introduced in the laser recrystallized silicon ribber by high-fluence 
low-energy (1.5 keV) ion implantation at a temperature of 250 °C. The deuterium con- 
centration profile presented in Fig. 1 was measured by secondary ion mass spectrometry 
(SIMS). As can be seen from Fig. 1, the analytical solution obtained provides good agree- 
ment with the experimental data [TT] . The following values of simulation parameters were 
used to fit the calculated curve to the experimental deuterium profile: Q = 4.32x10^^ 
cm-2; p^^ = 0.253; Rp = 0.035 //m; ARp = 0.028 /im; l^^ = 0.31 fim. Here Q is the 
fiuence of implanted hydrogen. Thus, it follows from the value of the fitting parameter p^^ 
that approximately 25% of the implanted deuterium atoms occupy interstitial positions. 
Migration of these nonequilibrium interstitial atoms results in the formation of extended 
"tail" on the deuterium concentration profile. 

In Fig. 2 the beryllium concentration profile after annealing calculated on the basis of 
analytical solution ([8]), (H7|) . (HHl) . and (H9l) is shown. For comparison, the experimental 
data of |12| are used. In Ref. [T^ beryllium diffusion in InGaAs/InP during rapid 
annealing at a temperature of 900 °C for 30 s was investigated. In the experiment under 
consideration, InP/InGaAs heterostructures were grown by gas-source molecular beam 
epitaxy onto semi-insulating (100) InP substrates. At first a 0.1 /.im InP buffer layer 
was grown, followed by 0.5 /im undoped InP and then 0.2 /im Be-doped Ino.53Gao.47As 
layer with a doping level of 3x10'' fim~^ was grown. Finally, an undoped InP layer of 
0.5 //m was grown on the Be doped layer. Then, a post-growth rapid thermal annealing 
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Figure 1: Deuterium concentration profile (solid line) calculated for the case of "hot" ion 
implantation at a temperature of 250 °C. The experimental data (dots) are taken from 
Sopori et al. [IT]. 

was performed in a halogen lamp furnace for 30 s. Beryllium profile measurements were 
made with SIMS. The measurements confirmed that the as-grown beryllium concentration 
profile is a step function, similar to the function (HTi) . 
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Figure 2: Beryllium concentration profile (solid line) calculated for the case of rapid 
thermal annealing (30 s at a temperature of 900 °C) of the InP/InGaAs heterostructures. 
The experimental data (dots) are taken from Ihaddadene-Lenglet et al. |12] . 
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To calculate the beryllium concentration profile after annealing it is supposed that 
during thermal treatment a part of substitutionally dissolved impurity atoms is transferred 
into interstitial positions. If the generation rate of beryllium interstitials is proportional 
to the impurity concentration, we can use function (HTj) to describe the distribution of the 
interstitial generation rate. It is also supposed that the impurity interstitial concentrations 
at the surface and in the bulk of a semiconductor are equal to zero. Then, the beryllium 
concentration profile after annealing is described by analytical expressions ([H), (HZj), (HHl) . 
and m^. 

As can be seen from Fig. 2, the calculated curve is in good agreement with the 
experimental data of [12]. The following values of simulation parameters were used to 
fit the calculated curve to the experimental concentration profile of beryllium: Cmax = 
2.672.x 10^ /im-3; p^^ = 0.254; xl = 0.5 /xm; xr = 0.7 /im; l^^ = 0.127 ^m. Here C^ax 
is the maximum concentration of substitutionally dissolved beryllium after annealing. 
It follows from the value of the fitting parameter that, as in the case of "hot" ion 
implantation, approximately 25% of the beryllium atoms occupy the transient interstitial 
positions. Migration of these nonequilibrium impurity interstitial atoms results in the 
formation of two "tails" on the beryllium concentration profile after annealing. 



5 Conclusions 

The analytical solutions of the equations describing impurity diffusion due to migration of 
nonequilibrium impurity interstitials are obtained. Two representative cases were inves- 
tigated: i) interstitial impurity diffusion during "hot" ion implantation under reflecting 
boundary condition at the surface of a semiconductor; ii) impurity interstitial diffusion 
from the doped epitaxial layer under the conditions of constant concentrations of inter- 
stitial impurity atoms at the surface and in the bulk. 

Using these solutions, we can verify the correctness of the approximate numerical 
calculations obtained by the codes intended for simulation of diffusion processes used in 
fabrication of semiconductor devices. Moreover, it is possible to carry out an analytical 
simulation of certain diffusion processes. As an example, hydrogen diffusion in silicon 
during high-fluence low-energy deuterium implantation at a temperature of 250 °C and 
beryllium diffusion from a doped epi-layer during rapid thermal annealing of InP/InGaAs 
heterostructures at a temperature of 900 °C were investigated. The impurity concentra- 
tion profiles calculated on the basis of the analytical solutions obtained agree well with 
the experimental data. Due to comparison with experiment, the values of the parameters 
describing interstitial impurity diffusion have been derived. It was obtained that for the 
processes under consideration approximately 25% of the impurity atoms occupied tran- 
sient interstitial positions. The average migration lengths of impurity interstitials are 0.31 
and 0.127 /xm for deuterium in silicon and beryllium in InP, respectively. 
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